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Abstract. Lithium and zirconium abundances (the latter taken as representative for s-process enrichment) are determined 
for a large sample of massive Galactic 0-rich AGB stars, for which high resolution optical spectroscopy has been obtained 
(R~40,000-50,000). This is done by computing synthetic spectra based on classical hydrostatic model atmospheres for cool 
stars using extensive line lists. The results obtained are discussed in the framework of "hot bottom burning" (HBB) and nu- 
cleosynthesis models. The complete sample is studied attending to various observational properties such as the position of the 
stars in the IRAS two-colour diagram ([12] - [25] vi- [25] - [60]), Galactic distribution, expansion velocity (derived from 
the OH maser emission) and period of variability (when available). We conclude that a considerable fraction of the sources 
observed are actually massive AGB stars (M > 3^ M©) experiencing HBB, as deduced from the strong Li overabundances 
found. A comparison of our results with similar studies carried out in the past for the Magellanic Clouds (MCs) reveals that, in 
contrast to MC AGB stars, our Galactic sample does not show any indication of s-process element enrichment. The differences 
observed are explained as a consequence of metallicity effects. Finally, we discuss the results obtained in the framework of 
stellar evolution by comparing our results with the data available in the literature for Galactic post-AGB stars and PNe. 
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1. Introduction 

The asymptotic giant branch (AGB) is formed by stars with ini- 
tial masses in the range between 0.8 and 8 Mq in a late stage of 
their evolution. The internal structure of an AGB star consists 
of an electron-degenerate C-O core surrounded by a shell of 
He. During most of the time H burning is the main source of 
energy for the AGB star but, occasionally, the inner He shell 
ignites in a "thermal pulse" and, eventually, the byproducts of 
He burning may reach the outer layers of the atmosphere in a 
"dredge-up" of processed material (the so-called third dredge- 

Send offprint requests to: D.A. Garcia-Hemandez 

* Based on observations obtained at the 4.2 m William Herschel 
Telescope operated on the island of La Palma by the Isaac Newton 
Group in the Spanish Observatorio del Roque de Los Muchachos of 
the Instituto de Astrofisica de Canarias. Also, based on observations 
obtained with the ESO 3.6 m telescope at La Silla Observatory (Chile) 



up). Repeated thermal pulses add carbon to the stellar surface. 
As a consequence, originally O-rich AGB stars can turn into 
C-rich AGB stars (C/O > 1). This scenario would explain the 
observed spectral sequence M-MS-S-SC-C in AGB stars (e.g. 
Mowlavi 1999). 

Another important characteristic of AGB stars is the pres- 
ence of neutron-rich elements (s-elements such as Sr, Y, Zr, 
Ba, La, Nd, Tc, etc.) in their atmospheres. These species are 
formed by slow-neutron captures in the intershell region. After 
a dredge-up episode, according to theoretical models, protons 
can be partially mixed in the '^C- and "^He-rich region between 
the hydrogen- and helium-burning shells (e.g. Straniero et al. 
1995, 1997; Mowlavi 2002; Lattanzio 2003), and can react with 
'^C to form '^C during the interpulse phase. The '^C(Q',n)'^0 
reaction releases neutrons, which are captured by iron nuclei 
and other heavy elements, forming s-elements that can later 
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be dredged up to the stellar surface in the next thermal pulse 
(Straniero et al. 1995, 1997; Busso et al. 2001). ^^Ne is an- 
other neutron source which needs to be considered. In this case 
the activation takes place during the convective thermal pulses 
(e.g. Straniero et al. 1995, 2000; Vaglio et al. 1999; GalUno 
et al. 2000). ^^Ne can be formed from '"^N and ^"*N is formed 
by the CNO cycle in the hydrogen shell. The ^^Ne(a,n)2^Mg 
neutron source requires higher temperatures and occurs typ- 
ically at higher neutron densities than the ^^C(Q',n)'*0 reac- 
tion. Thus, a different s-element pattern is expected depend- 
ing on the dominant neutron source. According to the most 
recent models the '^C(Q',n)^^0 reaction is the preferred neu- 
tron source for masses around 1-3 Mq, while for more massive 
stars (i.e. M > 3-4 Mq) neutrons are thought to be mainly re- 
leased through the ^^Ne(Q',n)^^Mg reaction (see, for example, 
Busso, GalUno & Wasserburg 1999 and Lattanzio & Lugaro 
2005 for a recent review). In the literature, there is strong ev- 
idence that most Galactic AGB stars enriched in s-process el- 
ements have masses around 1-3 Mq, where the ^^C(a,n)'^0 
reaction is the neutron donor (e.g. Lambert et al. 1995; Abia 
et al. 2001). Unfortunately, a confrontation of the predictions 
made by these models with observations of more massive AGB 
stars in our Galaxy is not yet available. 

In the case of the more massive O-rich AGB stars (M > 3^ 

Mq), the convective envelope can penetrate the H-burning shell 
activating the so-called "hot bottom burning" (hereafter, HBB). 
HBB takes place when the temperature at the base of the con- 
vective envelope is hot enough (7" > 2 X 10^ K) and '^C can be 
converted into ^^C and ^"^N through the CN cycle (Sackmann 
& Boothroyd 1992; Wood, Bessel & Fox 1983). HBB models 
(Sackmann & Boothroyd 1992; D'Antona & MazzitelH 1996; 
Mazzitelli, D'Antona & Ventura 1999) also predict the produc- 
tion of the short-lived ^Li by the chain ^He(Q',y)^Be (e",v)^Li, 
through the so-called "^Be transport mechanism" (Cameron & 
Fowler 1971). One of the predictions of these models is that 
lithium should be detectable, at least for some time, on the stel- 
lar surface. 

The activation of HBB in massive O-rich AGB stars is sup- 
ported by studies performed of AGB stars in the Magellanic 
Clouds (hereafter, MCs) (Plez, Smith & Lambert 1993; Smith 
& Lambert 1989, 1990a; Smith et al. 1995) which show a lack 
of high luminosity C stars beyond Mboi ~ -6. Instead, the more 
luminous AGB stars in the MCs are O-rich. The detection of 
strong Li overabundances together with strong s-element en- 
hancement in these luminous AGB stars in the LMC and in 
the SMC is the signature that these stars are indeed HBB AGB 
stars that have undergone a series of thermal pulses and dredge- 
up episodes in their recent past. The HBB nature of these stars 
is also confirmed by the detection of a very small '^C/'^C ra- 
tio (^ 3^), expected only when HBB is active (MazzitelU, 
D'Antona & Ventura 1999). Unfortunately, current HBB the- 
oretical models have been tested almost exclusively using the 
results obtained from the study of the more massive AGB stars 
in the MCs but have never been applied to Galactic sources, 
mainly because of the lack of observations available but also 
because of the inaccurate information on distances within our 
Galaxy. 



Indeed, only a handful of Li-rich stars have been found in 
our Galaxy so far (e.g. Abia et al. 1991, 1993; Boffin et al. 
1993) and, unlike those detected in the Magellanic Clouds, they 
are not so luminous (-6 < Mboi < -3.5). Most of them are 
low mass C-rich AGB stars (Abia & Isem 1996, 1997, 2000) 
and intermediate mass S- and SC-stars (Abia & Wallerstein 
1998) and not O-rich M-type stars. A few are less luminous 
red giant branch (RGB) stars at the bump (Charbonnel 2005). 
Some of these AGB stars are among the most Li-rich stars in 
our Galaxy, the so-called "super Li-rich" AGB stars (Abia et 
al. 1991). Under these conditions, however, HBB is not ex- 
pected to be active and the Li production is not well understood. 
Note that HBB is expected to be active only in the most mas- 
sive (and luminous) AGB stars (from ~4 to 7 Mq) (Mazzitelli, 
D'Antona & Ventura 1999), which should not be C-rich, but 
O-rich. Actually, the best candidates in our Galaxy are the 
so-called OH/IR stars, luminous O-rich AGB stars extremely 
bright in the infrared, showing a characteristic double-peaked 
OH maser emission at 1612 MHz. These stars are also known 
to be very long period variables (LPVs), sometimes with pe- 
riods of more than 500 days and large amplitudes of up to 2 
bolometric magnitudes. However, they experience very strong 
mass loss rates (up to several times 10"^ Mq yr"^) at the end 
of the AGB and most of them appear heavily obscured by thick 
circumstellar envelopes, making optical observations very dif- 
ficult (e.g. Kastner et al. 1993; Jimenez-Esteban et al. 2005a; 
2006a). Thus, no information yet exists on their lithium abun- 
dances and/or possible s-process element enrichment. 

In this paper we present the results obtained from a wide 
observational programme based on high resolution optical 
spectroscopy of a carefully selected sample of Galactic AGB 
stars thought from their observational properties to be massive. 
The criteria followed to select the sources in the sample are 
explained in Section 2. The high resolution spectroscopic ob- 
servations made in the optical and the data reduction process 
are described in Section 3, while the main results obtained are 
presented in Section 4. In this section we also show how spec- 
tral synthesis techniques were applied to derive the atmospheric 
parameters as well as the lithium and zirconium abundances of 
the stars in our sample. A discussion of these results in the con- 
text of HBB models, nucleosynthesis models and stellar evolu- 
tion can be found in Section 5. Finally, the conclusions derived 
from this work are given in Section 6. 

2. Selection of the sample 

A large sample of long period (300-1000 days), large ampU- 
tude variability (up to 8-10 magnitudes in the V band), late- 
type (> M5) O-rich AGB stars displaying OH maser emission 
with a wide range of expansion velocities (from just a few km 
s"' to more than 20 km s~') was carefully selected from the lit- 
erature. These sources are all very bright in the infrared and, as 
such, included in the IRAS Point Source Catalogue (Beichman 
et al. 1988). Some of them are so strongly reddened that they do 
not show any optical counterpart in the Digitized Sky Survey' 

' The Digitized Sky Survey was produced at the Space Telescope 
Science Institute under US Government grant NAG W-2166. 
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Fig. 2. IRAS two-colour diagram [12] - [25] vs [25] - [60], 
where the positions of all the sources in the sample are repre- 
sented (black dots). The continuous line is the "O-rich AGB 
sequence" (see text) which defines the sequence of colours ex- 
pected for O-rich AGB stars surrounded by envelopes with in- 
creasing thickness and/or mass loss rates (Bedijn 1987). 

(DSS) plates, although in a few cases the detection of these 
sources in the optical at a difl'erent epoch has been reported in 
the literature. 

Stars were included in the sample if they satisfied at least 
one of the above conditions (ideally as many of them as possi- 
ble), which guarantees that they are actually relatively massive 
stars. In Table 1 we list the 102 stars selected for the analysis, 
together with their IRAS names, other being names being taken 
from the literature, the galactic coordinates, the variability type 
taken from the Combined General Catalogue of Variable Stars 
(GCVS, Kholopov 1998), the spectral type taken from Kwok, 
Volk & Bidelman (1997) and references therein, and the IRAS 
colour indices [12] - [25] and [25] - [60]^. Similar informa- 
tion for a few well known Galactic M-supergiants and S-, SC- 
and C-type AGB stars that were also observed for comparison 
purposes is presented in Table 2. In addition, nine stars classi- 
fied as Mira-type stars in SIMBAD^ and included in the initial 
sample were found to exhibit observational properties corre- 
sponding to C-rich stars. These are listed in Table 3. 

The Galactic distribution of the stars in our sample is shown 
in Figure 1. The distribution observed clearly suggests that 
most of them must belong to the disc population. A few sources 
with bright optical counterparts located at high galactic lati- 
tudes may represent a subset of nearby stars. 

2 Defined as [12] - [25] = -2.5 logiFa/F.s) and [25]-[60] = -2.5 
log(F25/Fgo)' where F12, F25 and are the IRAS flux densities at 12, 
25 and 60 jum, respectively. 

^ Set of Identifications, Measurements, and Bibliography for 
Astronomical Data, operated at CDS, Strasbourg, France. 



In Figure 2 we show the position of these stars in the IRAS 
two-colour diagram [12] - [25] vs [25] - [60]. They fall, as 
expected and with very few exceptions, along the observational 
"O-rich AGB sequence" (Garcia-Lario 1992). This sequence 
agrees very well with the model predictions by Bedijn (1987) 
and is interpreted as a sequence of increasing thickness of the 
circumstellar envelopes and/or mass loss rates in O-rich AGB 
stars. 

3. Observations and data reduction 

The observations were carried out during several observing pe- 
riods in 1996-97. High resolution spectra were obtained us- 
ing the Utrecht Echelle Spectrograph (UES) installed at the 
Nasmyth focus of the 4.2 m William Herschel Telescope at the 
Spanish Observatorio del Roque de los Muchachos (La Palma, 
Spain) over three different observing runs in August 1996 (run 
#1), June 1997 (run #2) and August 1997 (run #3) and the 
CAsegrain Echelle SPECtrograph (CASPEC) of the ESO 3.6 
m telescope at the European Southern Observatory (La Silla, 
Chile) in February 1997 (run #4). The full log of the spectro- 
scopic observations is shown in Table 4, including more de- 
tailed information on the telescopes, dates of the observations, 
instrumentation used in each run, spectral dispersion, as well 
as on the spectral range covered. 

We used a TEK 1 124 x 1 124 CCD detector during the first 
run at the 4.2 m WHT using UES in August 1996 and a SITel 
2148 X 2148 CCD during the second and third runs in June 
1997 and August 1997, respectively. Since we were mainly in- 
terested in the spectral range between 6000 A and 8200 A, we 
used the 31.6 line/mm grating in order to provide full cover- 
age of this spectral region in a single spectrum. With a central 
wavelength around 6700 A our spectra extended over ~4000 
A during the first run and ~6000 A during the second and third 
runs. The spectra were taken with a separation between orders 
of around 21 pixels (or 7.5"), which is large enough to allow 
sky subtraction, taking into account that the targets in our sam- 
ple are non-extended. The resolving power was around 50 000, 
equivalent to a spectral resolution of 0.13 A around the Li I 
line at 6708 A. The selected set-up covered the spectral regions 
from 5300 A to 9400 A in about 47 orders with small gaps in 
the redder orders for the first run while the 4700-10300 A re- 
gion was covered in about 61 orders without any gaps for the 
second and third runs. 

CASPEC spectra, taken at the ESO 3.6 m telescope, cov- 
ered the wavelength range from 6000 A to 8200 A. The red 
cross-disperser (158 line/mm) was used, which gave a resolv- 
ing power of ~40000 (equivalent to a spectral resolution of 
~0.17 A around the Li I line at 6708 A) and an adequate in- 
terorder separation, using the TEK 1024 X 1024 CCD. With 
this set-up, the selected spectral region is completely covered 
in 27 orders with small gaps only between the redder ones. 

The observational strategy was similar in all runs. Since 
the stars of the sample are known to be strongly variable, tipi- 
cally 3-4 magnitudes in the /?-band and sometimes more than 
8 magnitudes in the V-band, we determined the exposure time 
according to the brightness of the source at the telescope. The 
typical exposure times ranged from 10 to 30 minutes. Orily for 




Fig. 1. Galactic distribution of the sources in the sample. 



a few very faint stars two 30 minute spectra were taken (and 
later co-added in order to increase the S/N ratio of the final 
spectrum). The goal was to achieve an S/N ratio of 50-150 in 
the region around 6708 A (the lithium line) in order to resolve 
this narrow absorption line, usually veiled at these wavelengths 
by the characteristic molecular bands which dominate the opti- 
cal spectrum of these extremely cool stars. Note that, because 
of the very red colours of the sources observed, the S/N ratios 
achieved for a given star can strongly vary from the blue to the 
red orders (e.g. 10-20 at -6000 A while >100 at -8000 A). 

At the telescope, all the 102 O-rich stars listed in Table 1 
were tried, but useful spectra were obtained only for 57 (56%). 
The remaining 45 (44%) sources were either too red to obtain 
any useful information on the strength of the Li 1 line at 6708 
A or the optical counterpart was simply not found, in the most 
extreme cases. Some targets with known, relatively bright opti- 
cal counterparts in the DSS plates were too faint to be detected 
on the date of the observations. In contrast, other targets with 
no optical counterpart in the DSS plates appeared as consider- 
ably bright stars. A few Galactic M-supergiants together with 
some C-, SC-, and S-type AGB stars were also observed for 
comparison (see list in Table 2). Some of the stars in this latter 
group are among the most Li-rich stars in our Galaxy. In addi- 
tion, as we have already mentioned in Section 2, a few stars in 
the initial sample turned out to be carbon stars when observed, 
and not O-rich as initially suspected (see list in Table 3). The 
total number of objects observed was 120. 

For wavelength calibration, several Th-Ar lamp exposures 
were taken every night during the UES runs. In the case of 
the CASPEC observations, these calibration lamp exposures 
were taken before or after any science exposure at the posi- 
tion of the target in order to keep control of possible dispersion 
changes with the telescope position. Note that with UES the 
system is more stable because the instrumentation is located at 
the Nasmyth focus of the 4.2 m WHT telescope. Finally, the 
corresponding bias and flatfield images were also taken at the 
beginning (or at the end) of the night. 

The two-dimensional frames containing the echelle spectra 
were reduced to single-order one-dimensional spectra using the 
standard echelle software package as implemented in IRAF.'* 

Image Reduction and Analysis Facility (IRAF) software is dis- 
tributed by the National Optical Astronomy Observatories, which 



Basically, the data reduction process consists of bias-level and 
scattered-light subtraction, the search and trace of apertures us- 
ing a reference bright star, the construction of a normalized flat- 
field image to remove pixel-to-pixel sensitivity fluctuations as 
well as the extraction of the 1-D spectra from the 2-D frames. 
For the wavelength caUbration we selected non- saturated Th- 
Ar emission lines and third or fourth order polynomials for the 
fitting. The calibration accuracy reached was always better than 
20 mA. Finally, we identified the terrestrial features (telluric 
absorption Unes) comparing the target spectra with the spec- 
trum of a hot, rapidly rotating star observed on the same night. 
It should be noted, however, that the majority of the spectral 
ranges used in the abundance analysis presented in this paper 
are not significantly affected by these features. 

4. Results 

4.1. Overview of the main spectroscopic properties 

In general, all stars show extremely red spectra with the flux 
level falling dramatically at wavelengths shorter than 6000 
A. In addition, the spectra are severely dominated by strong 
molecular bands mainly due to titanium oxide (TiO), as a con- 
sequence of the very low temperature and the O-rich nature of 
fliese stars.5 Thebandheads of TiO at -6651, 6681, 6714, 7055, 
and 7125 A are clearly present in all spectra. Interestingly, the 
bandheads of ZrO at -6378, 6412, 6474, 6495, 6505 and 6541 
A seem to be absent. These ZrO bandheads (together with those 
corresponding to other s-element oxides such as LaO or YO) 
are very strong in Galactic S-stars (see Figure 3). 

The TiO veiling effect is so intense that it is very difficult 
to identify individual atomic lines in the spectra of these stars 
with the exception of the Li I line at 6708 A, the Ca I lines 
at 6122 and 6573 A, tiie K I line at 7699 A, the Rb I Une at 
7800 A and a few strong Fe 1 lines. The K 1 and Rb 1 reso- 
nance fines sometimes show complex profiles, emission over 
absorption and blue-shifted components which may have their 

is operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

' In contrast, the spectra of the few C-rich AGB stars found in our 
spectroscopic survey (see list in Table 3) are dominated by CN and C2 
absorption bands. 
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IRAS name Othernames Galactic Coordinates Variability Spectral Type [12] -[25] [25] -[60] 
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origin in the expanding circumstellar shell. Actually, the dif- 
ference between the mean radial velocities of the stellar and 
circumstellar hue components is found to be of the order of the 
expansion velocity derived from the OH maser emission. In a 
few cases the Li I line can also show a similar behaviour (see 



Figure 4). Finally, some stars display Ha emission, which in 
this type of stars is generally interpreted as the consequence of 
the propagation of shock waves through the outer layers of the 
stellar atmosphere. 
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We detected the presence of the Li I resonance line at 6708 
A in 25% of the sources in the sample (25 stars) with a wide 
variety of strengths, while in 3 1% of them (32 stars) we did not 
find any signature of this line. The remaining 44% (45 stars) 
were either too red or the optical counterpart was simply not 
found at the moment of the observations. Sample spectra of 
stars showing increasing strength of the Li I line at 6708 A are 
presented in Figure 4. 



The observed sample has been divided into two different 
groups on the basis of the detection or non-detection of the 
lithium line in Tables 5 and 6, respectively. A third group is 
formed by the stars that were too red (or which did not show 
any optical counterpart at the telescope), and these are listed in 
Table 7. The IRAS name, the run in which the stellar field was 
observed, the OH expansion velocity and the pulsational period 
(if available) are given for every source in these tables. The 
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Table 2. The sample of comparison stars 



TT? AO? nri m - ^ 


Other names 


Galactic Coordinates 


Variability 


Spectral Type 


[12J-[Z3J 


tZDJ-[OUJ 


UUZio+joi / 


K Ana 


1 i / .U / 


T3 no 


Mira 


oj/4.je 


— U. / Zj 


—2. iU4 


m 1 no 1 1 


o r er 


1 J+.OZ 


1 1 Q 


CD ^ 


M4ela 


— U.4U / 


— i.sy / 




V Ant 

Zi Ant 


275.86 


+20.94 


oK 


OJ.4 


— U.ojU 


— Z.i /z 


18397+1738 


V821 Her 


47.78 


+10.01 


Mira 


Ce 


-0.873 


-1.499 


19133-1703 


TSgr 


20.32 


-13.01 


Mira 


S5/6e 


-1.125 


-1.372 


20166+3717 


WXCyg 


75.43 


+0.87 


Mira 


C9,2eJ 


-1.191 


+ 1.445 


23416+6130 


PZCas 


115.06 


-0.05 


SRc 


M3vla 


+0.070 


-1.538 


23554+5612 


WYCas 


115.57 


-5.62 


Mira 


S6/6e 


-0.633 


-1.286 


23587+6004 


WZCas 


116.76 


-1.91 


SRc 


SC7/10e 


-1.280 


-1.348 



Table 3. The sample of peculiar carbon AGB stars 



IRAS name 


Other names 


Galactic Coordinates 


Variability 


Spectral Type 


[12]-[25] 


[25]-[60] 


00247+6922 


V668 Gas 


120.86 


+6.87 


Mira 


C 


-0.774 


-1.762 


04130+3918 


RAFGL 6312 


160.41 


-8.09 




c 


-1.273 


-1.633 


09425-6040 


GLMP 260 


282.04 


-5.88 




c 


+0.790 


-1.045 


17297+1747 


V833 Her 


40.80 


+25.30 


Mira 


c 


-0.340 


-1.859 


19321+2757 


VI 965 Cyg 


62.57 


+3.96 


Mira 


c 


-0.703 


-1.598 


19594+4047 


V1968 Cyg 


76.52 


+5.59 


Mira 


c 


-0.285 


-1.594 


20072+3116 


VI 969 Cyg 


69.35 


-0.87 


Mira 


c 


-0.644 


-1.483 


23166+1655 


LLPeg 


93.53 


-40.35 


Mira 


c 


+0.101 


-1.238 


23320+4316 


LP And 


108.46 


-17.15 


Mira 


c 


-0.777 


-1.555 



Table 4. Log of the spectroscopic observations 



Run 


Telescope 


Instrumentation 


Dates 


Dispersion 

(A/pix) 


Spectral Range 

(A) 


#1 


4.2 m WHT 


UES 


16-17 August 1996 


0.065 


5300-9400 


#2 


4.2 m WHT 


UES 


10-12 June 1997 


0.065 


4700-10300 


#3 


4.2 m WHT 


UES 


18-20 August 1997 


0.065 


4700-10300 


#4 


3.6 m ESQ 


CASPEC 


21-24 February 1997 


0.085 


6000-8200 



same information is also presented for the comparison stars in 
Table 8 and for the sample of peculiar C-rich AGB stars in 
Table 9 (in this latter case the expansion velocities are derived 

from CO data). 

The detection of strong circumstellar Li I and Rb I lines 
in S Per (IRAS 02192+5821), sometimes classified in the lit- 
erature as an M-type Galactic supergiant, is remarkable. M- 
type supergiant stars are not expected to show a strong Li I 
line (e.g. Luck & Lambert 1982) or any s-process element en- 
hancement (e.g. Smith et al. 1995). According to our data, S 
Per looks like a possible genuine O-rich AGB star. In addi- 
tion, we find a cool effective temperature of ~3000 K which 
is quite different from Levesque et al. (2005) and further in- 
vestigation is needed. Sinoilarly, PZ Cas (IRAS 23416+6130), 
another source usually classified as an M-type supergiant (e.g. 
Levesque et al. 2005), is for the first time clearly identified 
by us as a possible S-,SC-type AGB star. PZ Cas seems to 
be slightly enriched in Zr and displays some strong Ba atomic 
lines but it is not enriched in Li, as expected if PZ Cas is a low 
mass S-,SC-type AGB star. Consistent with our possible AGB 
identification for the latter two sources, they show variability 



amphtudes of more than 3 magnitudes in the V-band and are 
the only two stars with periods well beyond 500 days in Table 
8. A more detailed analysis of these two stars is beyond the 
scope of this paper but it would be needed to reach a definitive 
conclusion. Also interesting is the new detection of a strong 
lithium line in the peculiar mixed chemistry (C-rich and O-rich) 
star IRAS 09425-6040, which is studied in a separate paper 
(Garcfa-Hemandez et al. 2006a). Lithium was also detected in 
the C-rich AGB stars IRAS 04130+3918, IRAS 20072+3116 
(V1969 Cyg) and IRAS 23320+4316 (LP And), for which a 
detailed analysis will be presented elsewhere. Finally, the S- 
star IRAS 10436-3459 (Z Ant), observed for the first time 
with high resolution spectroscopy, shows, as expected, intense 
molecular bands of s-element oxides (like ZrO and LaO) but no 
Uthium. This star wiU also be analysed in detail in a separate 
pubhcation. 

4.2. Modelling strategy 

Atmospheres of cool pulsating AGB stars present a major 
challenge for realistic, self-consistent modelhng. At present. 
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Fig. 3. High resolution optical spectra of sample stars displaying the lack of the ZrO absorption bands at 6474 and 6495 A com- 
pared with two Galactic S-stars (WY Cas and R And), where these bands are prominent. WX Ser (IRAS 15255+1944) and 
V697 Her (IRAS 16260+3454) are Li-detected while S CrB (IRAS 15193+3132) and IRAS 16037+4218 are Li undetected. The 
absorption band at ~6480 A corresponds to the TiO molecule. 



Pi 



IRASs 1 BU57 



\nAS 19129 



6700 

vr^i v^^ eri^l.Ti (A) 



Fig. 4. High resolution optical spectra of sample stars displaying increasing strength of the Li I line at 6708 A. Note the P 
cygni-type profile of the Li I line in IRAS 18057-2616. The jumps at 6681 and 6714 A correspond to bandheads of the TiO 
molecule. 



two types of self-consistent models exist. Classical hydrostatic Nordlund 1992). Time-dependent dynamical models for stellar 

model atmospheres include a sophisticated treatment of micro- winds (Winters et al. 2000; Hofner et al. 2003, and references 

physical processes and radiative transfer but neglect the effects therein) include pulsation, dust formation but a more crude de- 

of dynamics due to pulsation and winds (e.g. Plez, Brett & scription of radiative transfer. Propagating shock waves caused 
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Table 5. Galactic O-rich AGB stars where the lithium line was 
detected 







VexpiOH) 




Period 




IRAS name 


Run 


(km s-') 


Ref. 


(days) 


Ref. 


01085+3022 


3 


13.0 


3 


560 


1 


02095-2355 


1 










04404-7427* 


4 


7.7 


1 


534 


2 


05027-2158 


4 


7.9 


1 


368 


1 


05559+3825 


3 


single peak 


4 


590 


3 


06300+6058 


3 


12.3 


2 


440 


4 


07222-2005* 


4 


8.2 


1 


1200 


10 


07304-2032 


4 


7.4 


5 


509 


1 


09429-2148 


4 


12.0 


2 


650 


1 


11081-4203 


4 


single peak 


1 






11525-5057 


4 


single peak 


1 






12377-6102* 


4 


20.4 


1 






14337-6215 


4 


19.7 


9 






15211-4254 


4 


11.0 


1 






15255+1944 


1 


7.4 


2 


425 


1 


15576-1212 


4 


9.9 


2 


415 


1 


16030-5156 


4 


single peak 


1 






16260+3454 


1 


12.3 


2 


475 


1 


18025-2113 


1 


10.6 


1 






18057-2616 


1 


18.7 


1 






18413+1354 


2 


14.7 


4 


590 


3 


18429-1721 


1 


6.9 


1 






19129+2803 


2 


single peak 


1 


420 


10 


19361-1658 


2.3 


7.5 


1 






20052+0554 


2 


15.6 


12 


450 


3 


* Tentative Li detection 



References for the OH expansion velocities: 

1) te Lintel Hekkert et al. (1991); 

2) te Lintel Hekkert et al. (1989); 

3) Chengalur et al. (1993); 

4) Slootmaker, Habing & Herman (1985); 

5) Sivagnanam et al. (1989); 

6) Lewis, Eder & Terzian (1990); 

7) Lewis, David & Le Squeren (1995); 

8) Groenewegen et al. (2002) (from CO data); 

9) Sevensteretal. (1997); 

10) Groenewegen & de Jong (1998); 

11) Nymanetal. (1992); 

12) Lewis (1994) 

References for the periods: 

1) Combined General Catalogue of Variable Stars (GCVS), Kholopov 
(1998); 

2) Whitelock etal. (1994); 

3) Jones et al. (1990); 

4) Lockwood (1985); 

5) Engels et al. (1983); 

6) Nakashima et al. (2000); 

7) Slootmaker, Habing & Herman (1985); 

8) LeBertre(1993); 

9) Herman & Habing (1985); 

10) Jimenez-Esteban et al. (2006b) 
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Table 6. Galactic O-rich AGB stars where the lithium line was 
not detected 







V„,;(OH) 




Period 




IRAS name 


Run 


(km s-') 


Ref.i 


(days) 


Ref.i 


03507+1115 


3 


16.5 


2 


470 


1 


05098-6422 


4 


6.0 


2 


394 


2 


05151+6312 


1 


14.5 


2 




3 


07080-5948 


4 


5.0 


1 






07445-2613 


4 


4.8 


2 


391 


1 


10189-3432 


4 


3.4 


2 


303 


1 


10261-5055 


4 


3.7 


1 


317 


1 


13379-5426 


4 


10.2 


1 






13442-6109** 


4 


23.7 


1 






13475-4531 


4 


3.7 


1 






14086-0730* 


4 


13.7 


2 






14086-6907* 


4 


13.4 


1 






14247+0454 


4 


4.0 


2 


354 


1 


14266-4211 


4 


9.1 


1 






15193+3132 


1 


3.3 


2 


360 


1 


15586-3838 


4 


9.3 


1 






16037+4218 


1 


3.5 


1 


360 


10 


16503+0529 


2 


3.8 


3 


323 


1 


17034-1024 


2 


single peak 


1 


346 


1 


17359-2138* 


2 


8.5 


1 






18050-2213 


3 


19.8 


2 


732 


1 


18304-0728 


3 


15.3 


2 






18454-1226** 


3 


single peak 


1 






19147+5004 


1 


12.6 


1 






19157-1706 


2 


14.7 


1 






19412+0337 


2 


14.5 


6 


440 


3 


20343-3020 


3 


8.4 


1 


349 


2 


20350+3741 


3 






450 


3 


22180+3225 


1 


4.3 


3 






22466+6942 


1 


5.6 


1 


408 


6 


22510+3614 


1 


2.8 


3 






23492+0846** 


1 











* Poor S/N ratio 

** Possibly not an AGB star 



* See applicable references in Table 5 



by stellar pulsation modify the structure of the atmosphere on 
local and global scales, which give rise to strong deviations 
from a hydrostatic stratification. The radiative field is domi- 
nated by the effect of molecular opacities or even by the dust 
grains forming in the cool outer layers of these atmospheres. 
Chemistry and dust formation may be severely out of equihb- 
rium. 

Hydrostatic model atmospheres based on the MARCS code 
(Plez, Brett & Nordlund 1992; Gustafsson et al. 2003) repro- 
duce well the observed spectra of C-rich AGB stars (e.g. Loidl, 
Lancon & Jcergensen 2001) and O-rich AGB stars (e.g. Alvarez 
et al. 2000) in the optical domain and were thus adopted for the 
analysis. Indeed, the optical spectra of massive O-rich AGB 
stars in the Magellanic Clouds have succesfully been modelled 
using this code (e.g. Plez, Smith & Lambert 1993). A compar- 
ison of these models with high resolution observational data 
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Table 7. Galactic O-rich AGB stars for which no analysis of 
the Uthium line was possible 



Table 8. Comparison stars 



„„(OH) 
') 



TR A S namp 

1 1 V /\ wj> 11 till It 


Run 




nin'?7-i-i9io 


1 


1 8 9 






1 1 


\tZ.J 1 Ut^Ut^J J 




1 S S 

U.J 






1 6 


(16907-1- ADAS 




1 9 


(10104— 4"^ 18 


4 


1 1 ^ 

1 1 .J 




4 


1^9 

1 J.Z 


1 ^9(1^ — S^^fi 


4 


1 9 7 
iz. / 




4 


9S 1 

ZJ. 1 


1 ^"^4.1 —694.fi 


4 


16 7 


1 Jji / OJiJ 


4 


17 
1 1 .y 




4 


1 8 8 
1 o.o 




4 


94 


1 7 1 (1^— (l^SQ 


9 ^ 


1 J. J 


179^0—9819 




90 7 


174-1'J_17<(1 


-J 

J 


1 


174^^-9^9^ 




1 
ly.y 


1744"^- 9S10 




16 1 

lU. i 


i7';ni — 9fi'^fi 

1 / J\J I — Z.UJU 


-1 


9"? 


18(171—1797 


T. 

J 


1 1 


1 8083-9630 


9 3 


18.2 


1 8 1 79_9^n^ 




1 S 


18108-1940 




1 5 ^ 


189'^7— 1000 


9 


10 1 

17.1 


1 897"?- 07'?8 


9 

z 


14 ^ 


1 8976-1431 


2 


12.0 


1 8^1 9-1 900 




19 3 

IZ. J 


1 8^14—1 1^1* 






1 8'?48— 0'>96 


9 

Z, J 


1"? 


1 84^9—0140 


9 ^ 

Z, J 


17 


1 84^7—064^ 


9 ^ 

Z, J 


19 1 

IZ. i 


1 8460— 09 S4 


9 ^ 

Z, J 


90 9 
zu.z 


18488—0107 
lotoo WIU/ 


-J 


90 'i 

ZW. J 


1 8S40-I-0708 




1^ 7 


1 8'^60-^06^8 


1 9 ^ 

i ,Z, J 


1 6 


inncQ 00 1Q 


9 ^ 

zp 


1 J. J 


19161+2343 


2,3 


17.9 


19192+0922 


2,3 


16.5 


19254+1631 


2 


19.0 


19426+4342 


1,3 


8.6 


20077-0625 


1,3 


12.3 


20109+3205 


3 


5.7 


20181+2234 


1,3 


11.2 


20272+3535 


3 


12.0 


22177+5936 


3 


15.1 



Ref.' 



Period 

(days) Ref.' Notes 



420 



928 
1488 



10 



2 


660 


2 


R 


2 


1994 


7 




7 


534 


3 


R 


2 


635 


1 


R 


2 


520 


10 


R 








R 








R 



R 
R 



R 
R 

R 



845 


9 




1975 


9 




700 


5 


R 


890 


9 


R 






R 


1575 


3 




1140 


5 




652 


9 




1730 


5 




1540 


4 




840 


9 




1033 


5 




510 


1 


R 


552 


9 




1162 


9 


R 






R 


680 


1 


R 


382 


1 




1603 


9 




1215 


9 





* Single peak OH maser 

R: optical counterpart too red; 

. . .: optical counterpart not found 

' See applicable references in Table 5 



in the optical for massive O-rich AGB stars in our Galaxy is 
unfortunately still lacking. 







v„„(OH) 




Period 


IRAS name 


Run 


(kms-') 


Ref.i 


(days) Ref.' 


00213+3817 


1 


9.0 


10 


409 1 


02192+5821 


1 


14.5 


2 


822 1 


10436-3459 


4 


11.3 


1 


104 1 


18397+1738 


2 


13.4 


8 


511 3 


19133-1703 


1 


10 


10 


395 1 


20166+3717 


3 






410 1 


23416+6130 


1 


26.5 


2 


925 1 


23554+5612 


1 


15.5 


10 


477 1 


23587+6004 


1 






186 1 



' See applicable references in Table 5 



Table 9. Peculiar C-rich AGB stars 



IRAS name 


Run 


v„„(CO) 
(km s"') 


Ref.' 


Period 
(days) 


Ref.' 


Notes 


00247+6922 


3 


17.2 


8 


650 


3 




04130+3918 


3 






470 


10 


B 


09425-6040 


4 










B 


17297+1747 


2 


15.2 


11 


520 


3 


R 


19321+2757 


2 


24.4 


8 


625 


3 


R 


19594+4047 


3 


20.5 


8 


783 


3 




20072+3116 


2 


25.6 


8 


550 


3 


R 


23166+1655 


3 


15.1 


8 


520 


3 




23320+4316 


3 


14.7 


8 


620 


3 


R 



B: bright optical counterpart; 
R: optical counterpart too red; 
. . .: optical counterpart not found 

' See applicable references in Table 5 



Our analysis combines state-of-the-art line-blanketed 
model atmospheres and synthetic spectroscopy with extensive 
Une Usts. For this we need to estimate first the range of values 
of the stellar parameters that can reasonably be adopted for the 
stars in our sample: effective temperature T^g, surface gravity 
log g, mass M, metallicity Z = [Fe/H], microturbulent velocity 
^ and C/O ratio. 

Spherically synraietric, LTE, hydrostatic model atmo- 
spheres were calculated using the MARCS code for cool stars 
_(Gustafsson et al. 2003). The models are designed with the no- 
tation (Jeff, log g, M, Z, C/O), where Jeff is the efi'ective 
temperature of the star in K, g is the surface gravity in cm s"^, 
M is the stellar mass in M©, Z is the metallicity, f is the micro- 
turbulent velocity in km s"' and C/O is the ratio between the C 
and O abundances. 

Synthetic spectra were generated with the 
TURBOSPECTRUM package (Alvarez & Plez 1998) which 
shares much of its input data and routines with MARCS. 
Solar abundances of Grevesse & Sauval (1998) were always 
adopted, except for the iron abundance, which was taken as log 
e(Fe)= 7.50. Most opacity sources for cool M-type stars were 
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included and taken from the literature. In addition, updated 
line lists were produced by us for TiO (Plez 1998), ZrO (Plez 
et al. 2003) and VO (Alvarez & Plez 1998). 

For atomic lines, the primary source of information was the 
VALD-2 database (Kupka et al. 1999). The NIST'' atomic spec- 
tra database was also consulted for comparison. Where pos- 
sible, gf values were checked by fitting the solar spectrum. 
VALD-2 gf values taken from Kurucz (1993, 1994) were often 
erroneous and did not fit the atomic lines of the solar spectrum. 
The identification of features was made using the solar atlases 
by Moore, Minnaert, & Hootgast (1966) and Wallace, Hinkle 
& Livingston (1993, 1998) and the observed solar spectrum 
of Neckel (1999). The TURBOSPECTRUM program was run 
using the abundances from Grevesse & Sauval (1998) together 
with the solar model atmosphere by Holweger & Miiller ( 1 974) 
with parameters T^s - 5780 K, log g - 4.44 and a variable mi- 
croturbulence as a function of the optical depth. 

The whole machinery was tested on the high resolution op- 
tical spectrum (R ~ 150000) of the K2 IIIp giant Arcturus 
(a Boo) from Hinkle et al. (2000). A MARCS model atmo- 
sphere was used with the fundamental parameters determined 
by Decin et al. (2003a): T^s = 4300 K, log g = 1.50, Z = -0.5, 
M = 0.75 Mo, ^ = 1.7 km s"', C/N/O = 7.96/7.55/8.67 and 
'^C/'^C = 7. These parameters are in excellent agreement with 
other determinations reported in the literature. Some Kurucz 
gf values of metallic lines of Fe I, Ni I, Co I, Cr I, Ti I, Zr I 
and Nd II, which were too weak in the solar spectrum, were 
adjusted so as to yield a good fit to the Arcturus spectrum in 
several regions (especially in the 7400-7600 A and 8 1 00-8 150 
A regions); otherwise we used the gf values from the VALD-2 
database. This exercise was very useful in confirming the lack 
of s-process atomic lines in the 7400-7600 and 8100-8150 
A spectral windows in the Galactic O-rich AGE stars in our 
sample, as we shall see later 

4.3. Spectral regions of interest 

For the abundance analysis we concentrated our attention on 
the following spectral regions: 

4.3.1. Lithium 

We used the Li I resonance line at -6708 A to derive the lithium 
abundances, which can be used as a signature of HBB. The se- 
lected spectral region ranges from 6670 to 6730 A and cov- 
ers the TiO molecular bandheads at ~6681 and 6714 A, which 
are sensitive to variations in the effective temperature. In ad- 
dition, the spectral regions around the subordinate Li I lines 
at ~6103 and 8126 A were also inspected. However, we de- 
cided to drop them from the analysis because the latter lines 
are much weaker Moreover, the S/N ratio is usually very low 
around 6103 A in our sources and the contamination by telluric 
lines at 8126 A is very strong. 



* National Institute of Standards and Technology (NIST); 
|http://physics.nist.gov/cgi-bin/AtData/lines_form, 



4.3.2. Rubidium and Potassium 

We also synthesized the spectral regions around the resonance 
lines of Rb I at -7800 A and K I at 7699 A, with the inten- 
tion of using the elemental abundances derived from these lines 
as neutron density and metallicity indicators, respectively. In 
addition, these spectral regions were also used to check that 
the model parameters adopted to reproduce the lithium region 
also provided a reasonably good fit at other wavelengths. Two 
intervals covering ~60 A in the regions 7775-7835 A and 
7670-7730 A were selected for this purpose. Unfortunately, 
the frequent detection of circumstellar components to these 
lines prevented the accurate determination of the K I elemental 
abundances, and they will not be discussed here. To a lesser ex- 
tent the problem also affects the Rb I line. A detailed analysis 
of the Rb abundances will be presented in a forthcoming paper 
(Garcia-Hernandez et al. 2006b). 

4.3.3. s-process elements 

The study of the s-process elements was carried out through 
the analysis of several ZrO absorption bands in the region 
from 6455 to 6499 A. In particular, we used the ZrO band- 
head around ~6474 A to determine severe upper limits to the 
Zr elemental abundance, taken as representative of the overall 
s-process element enrichment. Surprisingly, this band was not 
detected in any star in our sample, neither did we find any sig- 
nature from atomic Zr or from atomic lines corresponding to 
any other abundant s-process element in the 7400-7600 A and 
8100-8150 A spectral windows, where many atomic lines of 
Zr I, Nd II, Ba II might potentially have been found, as we 
show later. 

4.4. Derivation of thie stellar parameters 

In order to analyze how the variations in stellar parameters in- 
fluence the output synthetic spectra, we constructed a grid of 
MARCS model atmospheres and generated the associated syn- 
thetic spectra with the following specifications: i) the stellar 
mass was in all cases taken to be 2 M©; ii) Teff values ranging 
from 2500 to 3800 K in steps of 100 K; iii) C/O ratio values 
of 0.5, 0.7, 0.8 and 0.9; iv) log g between -0.5 and 1.6 dex in 
steps of 0.3 dex; v) microturbulent velocity ^ between 1 and 
6 km s ' in steps of 0.5 km s"'; vi) metallicity Z between 0.0 
and -0.3 dex; vii) log e(Zr)^ from 1.6 to 3.6 dex in steps of 
0.25 dex; viii) CNO abundances shifted +1.0 dex (in steps of 
0.5 dex) from the solar values of Grevese & Sauval (1998); and 
ix) '^C/'^C ratios of 10 (as expected for HBB stars) and 90 (the 
solar value). Finally, the synthetic spectra were convolved with 
a Gaussian profile (with a certain FWHM typically between 
200 and 600 mA) to account for macroturbulence as well as 
instrumental profile effects. 

This, however, results in an extremely large grid of syn- 
thetic spectra containing thousands of possible combinations of 
the above parameters ! In order to reduce the number of spectra 

^ The zirconium abundance is given in the usual scale 12 + log 
A'(Zr). On this scale the solar zirconium abundance is 2.6 dex. 
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to be considered in our analysis, further constraints were im- 
posed by studying the sensitivity of our spectral synthesis to 
changes in these parameters. 

4.4.1 . Stellar mass 

As we have already mentioned, the stellar mass was in all cases 
considered to be constant and equal to 2 Mq . This is because the 
temperature and pressure structure of the model atmosphere is 
practically identical for a 1 Mq and a 10 Mq model atmosphere 
(see figure 1 in Plez 1990). Increasing the mass wiU just have 
a marginal effect on the output synthetic spectrum through a 
decrease, for a given gravity, of the atmospheric extension. 

4.4.2. Effective temperature 

The synthetic spectra are, in contrast, particularly sensitive to 
Teff , which determines the overall strength of the molecular ab- 
sorption over the continuum (mainly, from TiO molecules but 
also from VO, and from ZrO if the Zr elemental abundance is 
increased above a certain Umit). Thus, we decided to keep the 
whole range of values initially considered in our grid, i.e. Teff 
from 2500 to 3800 K in steps of 100 K. Actually, the depth 
of the molecular bands increases considerably with decreasing 
temperature. The effective temperature has also a large impact 
on the strength of the Li I line at 6708 A in the synthetic spec- 
trum. A decrease of Teg has the effect of increasing the TiO 
veiling, as well as the strength of the Li 1 line, as a consequence 
of the changes in the Li 1/Li 11 population equilibrium. 

4.4.3. C/0 ratio 

The C/O ratio basically determines the prevalence of O-rich 
molecules C/O < 1) against C-rich molecules like CN, C2, etc. 
In the case of our O-rich AGB stars, the precise value of the 
C/O ratio adopted, always < 1, affects mainly the strength of 
the TiO veiling and of the Li 1 line in the synthetic spectrum. 
An increase of the C/O ratio will decrease the TiO veiling and 
increase the Li 1 line strength relative to the adjacent contin- 
uum in a very similar way as a decrease in Tgg. This is a good 
example of how a set of parameters providing a good match to 
the observations for a given spectral region is not necessarily 
unique. A good fit may be obtained for another combination 
of parameters that can correspond to quite different Li abun- 
dances. Fortunately, one of these two sets of parameters usu- 
ally does not provide acceptable fits for other spectral regions. 
In particular, the presence of detectable atomic fines in the syn- 
thetic spectrum is very sensitive to variations in the C/O ratio: 
an increase in the C/O ratio will lower the TiO veiling, making 
the detection of atomic fines easier. According to the models, 
the immediate effect should be the detection of much stronger 
atomic lines of K I, Fe I, Zr I, Nd 11, etc. Since we do not see 
any such effect in our spectra, we conclude that C/O must al- 
ways be < 0.75 in our stars. Abundances derived from models 
with C/O between 0. 15 and 0.75 show actually little differences 
(Plez, Smith & Lambert 1993). Given that all stars in our sam- 
ple are clearly O-rich and taking into account the above consid- 



erations, we decided to use C/O - 0.5 in all cases, as a constant 
value. This selection is in agreement with other more detailed 
determinations made in massive O-rich AGB stars in the MCs 
(e.g. Plez, Smith & Lambert 1993; Smith et al. 1995) and in 
a few low mass M-type stars studied in our Galaxy (Smith & 
Lambert 1985, 1990b).^ 



4.4.4. Surface gravity 

The surface gravity also affects the appearance of the output 
synthetic spectra but its effect is small compared to that of the 
effective temperature. The influence of a change in the value 
adopted for the surface gravity within the range of values (log 
g between -0.5 and 0.5 dex) under consideration is not as se- 
vere as a temperature change. The molecular absorption be- 
comes slightly weaker at higher surface gravities. For example, 
an increase in the surface gravity of 0.5 dex has approximately 
the same impact as an increase of 100 K in effective tempera- 
ture. However, for a fixed temperature we found that the lowest 
gravities in general fit better both the TiO band strength and 
the pseudo-continuum around the 6708 A Li I fine. As the ap- 
pearance of the spectra is not so dependent on the surface grav- 
ity, and considering that its value must be low in these mass- 
losing stars, a constant surface gravity of log g = -0.5 was 
adopted for all the stars in the sample. This selection seems to 
be also appropiate if it is compared with the values adopted by 
Plez, Smith & Lambert (1993) of -0.33 < log g < 0.0 for O- 
rich AGB stars in the MCs with effective temperatures between 
3300 and 3650 K. 



4.4.5. Microturbulent velocity 

The microturbulent velocity is usually derived by demanding 
no correlation between the abundance of iron derived from in- 
dividual lines and their reduced equivalent widths. Previous 
studies of O-rich AGB stars found microturbulent velocities be- 
tween 2 and 4 km s"^ (Plez, Smith & Lambert 1993; Smith & 
Lambert 1985, 1989; Smith et al. 1995; Vanture & Wallerstein 
2002). Unfortunately, we cannot estimate the microturbulent 
velocity in the stars of our sample due to the lack of useful Fe 
1 atomic lines (the small number of detectable Fe I lines does 
not cover a wide range of equivalent widths). Thus, we will as- 
sume in the following a microturbulent velocity f = 3 km s~', 
which is a typical value generally adopted for AGB stars (e.g. 
Aringer, Kerschbaum & J0rgensen 2002). Higher values of ^ 
would strenghten all atomic lines in the synthetic spectra and 
we do not see this effect in our data. 



* This is the only other chemical analysis previously carried out in 
Galactic O-rich AGB M-type stars (only six stars!) before the work 
presented here. However, these were low mass stars (M < 1.5-2Mq) 
and their study was focussed on the determination of CNO abundances 
through near-IR spectra and the s-process enrichment, but they did not 
study the lithium enrichment in particular. 
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4.4.6. Metallicity 

The metallicity of the Galactic O-rich AGB stars in our sample 
is assumed to be solar, as expected for stars belonging to the 
disc population of our Galaxy. A lower metallicity is unlikely 
since all the stars are expected to be at least of intermediate 
mass (M > 2-3 Mq), if not more massive. This assumption is 
also in good agreement with the typical metallicities derived for 
other AGB stars in our Galaxy (e.g. Abia & Wallerstein 1998; 
Vanture & Wallerstein 2002). The Ca I Unes at -6122 A and 
6573 A generally used for this determination in other studies 
are unfortunately not sensitive enough to metallicity nor to sur- 
face gravity at the very low temperature of our sample stars (see 
figure 7 of Cenarro et al. 2002). This was checked by running 
different test models and spectral syntheses in the 6100-6160 
A and 6535-6585 A regions where these lines fall. A decrease 
in the metallicity implies less availability of metals like Ti, and 
thus, lowers the TiO veiling, making the detection of atomic 
lines easier. Actually, this may be the main reason why they 
are easily detected in O-rich AGB stars at the metallicity of the 
MCs (e.g. Z = -0.5 in the SMC, Plez, Smith & Lambert 1993), 
while we do not detect them in our spectra. 

4.4.7. Zirconium 

The synthetic spectra are also very sensitive to changes in 
the zirconium elemental abundance, especially in those wave- 
length regions where ZrO molecular bands are present. For a 
set of fixed stellar parameters, an increase in the zirconium 
abundance leads to stronger bands. Since the zirconium abun- 
dance was a priori unknown in our stars, we decided to use a 
large set of zirconium abundances with log e(Zr) between 1 .6 
and 3.6 dex in steps of 0.25 dex. The efi'ect of a variation in the 
zirconium elemental abundance on the ZrO bands is in general 
much stronger than any other efTect associated with variations 
in the T^g or C/O ratio. In contrast, even at relatively high zir- 
conium abundance, the corresponding Zr atomic lines are more 
diflicult to detect at low T^ff or low C/O ratio because the TiO 
and ZrO lines become stronger and more numerous under these 
conditions compared to the Zr atomic lines. 

4.4.8. CNO abundances and ^^CPC ratio 

In the synthetic spectra, relative variations in the carbon, ni- 
trogen and oxygen abundances and of the '^C/'^C ratio are 
completely negligible with respect to other stellar fundamen- 
tal parameters in the spectral regions under analysis. Thus, we 
decided to fix them and adopt their solar values in the following 
analysis. 

4.5. Abundance determination 
4.5.1. Overall strategy 

After imposing these further constraints, the grid of MARCS 
model spectra to consider is composed of "only" a few hun- 
dred spectra, with efi'ective temperatures in the range T^g = 
2500-3800 K in steps of 100 K, log e(Zr) between 1.6 and 3.6 



dex in steps of 0.25 dex, and a variable value of the FWHM 
in the range 200-600 mA in steps of 50 mA, keeping all the 
other stellar parameters fixed: log g = -0.5, M = 2 Mq, solar 
metallicity Z =[Fe/H]=0.0, ^ = 3 km s ', C/O = 0.5, solar CNO 
abundances and solar '^C/^^C ratios. 

In order to find the synthetic spectrum which better fits the 
observed spectrum of a given star a modified version of the 
standard test was used. When fitting observed data Yobsi 
to model data Ysynth-,, the quality of the fit can be quantified 
by the test. The best fit corresponds to that leading to the 
minimum value oix^- This test is used here in a modified way 
and is mathematically expressed as: 

N 2 

2 _ Y' [yobsj - Ysynthi{x\_...XM)\ 
^ ~ Zj Yobsi ^ ' 

with the number of data points, and M the number of free 
parameters. 

The observed spectra, once fully reduced, were shifted to 
rest-wavelength using the mean radial velocity shift derived 
from the Li I and Ca I Unes. In addition, they were re-binned to 
the same resolution as the synthetic ones (0.06 A/pix) and nor- 
maUzed in order to make the comparison easier. The observed 
spectra were then compared to the synthetic ones in intervals 
of ~60 A which allowed the analysis of their overall character- 
istics as well as of the relative strength of the TiO molecular 
bands. 

The fitting procedure makes a special emphasis on the 
goodness of the fit in the hthium region (6670-6730 A). In this 
spectral range the goal was to fit the relative strength of the TiO 
bandheads at ~6681 and 6714 A, which are very sensitive to 
the eff'ective temperature, together with the pseudo-continuum 
around the Li I line. We first determined by x^ minimization 
which of the spectra from our grid of models provided the best 
fit, mainly fixing T^s. The best fits resulting from the;^''^ min- 
imization were also judged by eye in order to test the method. 
The lithium content was then estimated by changing the lithium 
abundance. This procedure was repeated on each star of the 
sample for which an acceptable S/N ratio was achieved around 
the Li 1 6708 A line. Unfortunately, we could not analyse a few 
AGB stars with a low signal-to-noise ratio in their spectra at 
6708 A (see Table 7). IRAS 18025-2113, IRAS 03507+1115 
and IRAS 18304-0728 showed unusual spectra, very differ- 
ent from the rest of the stars in our sample and a good fit 
from our grid of MARCS model spectra was not possible. The 
spectra of these stars show unusual TiO band strengths which 
could only be reproduced with a low T^s but the rest of the 
spectrum (e.g. the local continuum level) seems to be hotter. 
In addition, the atomic lines are broader (IRAS 18025-2113 
and IRAS 18025-21 13) or narrower (IRAS 03507+1 1 15) than 
the FWHM needed to describe the TiO bandheads and Unes. 
Their complicated spectra look like a combination of two tem- 
peratures and FWHMs suggesting that they could probably 
be double-lined spectroscopic binaries or strongly affected by 
shock waves propagating in their atmosphere. 

The best model fit in the Li I spectral region (6670-6730 
A) is also usually found to fit the wavelength regions around 
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Fig. 5. Best model fit and observed spectrum around the Li I line (6708 A) for the star IRAS 11081-4203. The Li abundance 
derived from this spectrum was log e(Li) =1.3 dex. The synthetic spectra obtained for Li abundances shifted +0.2 dex and -0.2 
dex from the adopted value are also shown. The parameters of the best model atmosphere fit are indicated in the top label. 
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Fig. 6. Best model fit and observed spectrum in the region 6670-6730 A for the star IRAS 11 08 1 -4203. The T^ff derived from 
this spectrum was 3000 K. The parameters of the best model atmosphere fit are indicated in the top label. 



the ZrO bandhead and the K 1 and Rb 1 resonance lines rea- 
sonably well. In general, the overall shape of the spectrum (in- 
cluding the TiO bandheads) is very well reproduced. As an ex- 
ample, the fits made in diff'erent spectral regions for the star 
IRAS 11081-4203 are presented in Figures 5 to 9. The eff^ec- 
tive temperatures of the best-fitting model spectra together with 



the lithium abundances (or upper limits) derived following the 
above procedure are listed in Table 10, where we have sepa- 
rated the Li detected stars from the Li undetected ones. 

For the measurement of the zirconium abundance we con- 
centrated our attention on the synthesis of the 6455-6499 
A spectral region. In this case, the goal was to obtain the best 
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Fig. 7. Best model fit and observed spectrum in the region 6460-6499 A for the star IRAS 11 08 1 -4203. The Zr abundance derived 
from this spectrum was log e(Zr) = 2.6 dex, which corresponds to the solar value. The parameters of the best model atmosphere 
fit are indicated in the top label. Note the non-detection of the characteristic strong ZrO bandheads typical of Galactic S-stars. 



fit around the ZrO molecular bands at 6474 A and 6495 A. The 
non-detection of these features in any of the stars under analy- 
sis imposes severe upper limits to the zirconium abundance. To 
be sure that there was no other way to interpret our data, we also 
checked the presence of atomic zirconium in the 7400-7600 
and 8100-8150 A spectral regions, where strong Zr 1 atomic 
lines should have been detected in that case. Unfortunately, the 
local continuum in these latter wavelength regions is not well 
reproduced in the O-rich AGB stars in our sample, in contrast 
with the perfect agreement obtained in some of the compari- 
son stars observed, such as the S-star IRAS 10436-3454, for 
which a perfect fit is obtained both for the ZrO bandheads at 
6455-6499 A (see Figure 12) as well as for the Zr I atomic 
lines in the 7400-7600 and 8100-8150 A spectral regions (not 
shown). Sample spectra of the region around 7565 A are shown 
in Figure 10, where the position of some atomic lines of Zr I, 
Nd II and Fe I are indicated. As we can see, the atomic lines 
corresponding to s-process elements (e.g. Zr, Nd, La, etc.) were 
not detected in any of the sample stars observed. In contrast, 
these s-element atomic lines appear very strong in the Galactic 
S-star IRAS 10436-3454. Note, however, that the stars in our 
sample are the coolest O-rich AGB stars ever studied (with ef- 
fective temperatures between 2700 K and 3300 K; see Table 
10). We suspect that the strong discrepancies observed between 
model and observations may be due to the non-inclusion of 
other O-rich molecules, such as H2O, in our line lists. The ef- 
fect of H2O may be ignored in S-stars, where the C/O ratio 
is close to unity, but it might become dominant in our O-rich 
stars at these very low temperatures, preferentially at the longer 
wavelengths (e.g. Allard, Hauschildt & Schwenke 2000; Decin 
et al. 2003a, 2003b). In addition, lanthanum oxide also has a 



strong absorption band around 7500 A and it could have been a 
good candidate for detection in the spectral range here consid- 
ered. As for ZrO, the non-detection of the strong LaO molec- 
ular bandhead, usually very strong in S-stars at these wave- 
lengths, is another indication of the lack of s-process elements 
in our sample stars. 



4.5.2. Lithium and s-element abundances 

The derived Li abundances are displayed in Table 10, where 
we can see that almost all Li detected stars show enhanced Li 
abundances log e(Li) > 1, i.e. larger than solar, but smaller than 
those found in the so-called "super Li-rich" AGB stars (with 
log e(Li) > 3-4; e.g. Abia et al. 1991) in our Galaxy. A very 
similar range of Li overabundances was found in the massive 
O-rich AGB stars studied in the MCs (Plez, Smith & Lambert 
1993; Smith & Lambert 1989, 1990a; Smith et al. 1995). The 
errors in the derived Li abundances mainly reflect the sensitiv- 
ity of our models to changes in the adopted stellar parameters. 
In particular, they strongly depend on the uncertainties in the 
determination of the effective temperature (+100-200 K; larger 
for the coolest stars) and metallicity (±0.3), while they are less 
sensitive to the microturbulent velocity (±1 km s"'), surface 
gravity (±0.5) and FWHM (±50 mA) uncertainties. The ef- 
fect of the uncertainty in the location of the pseudo-continuum 
around the Li I line on the measured abundances is negligi- 
ble compared with any other change in the adopted stellar pa- 
rameters. If we consider all these uncertainties as independent 
sources of error, the resulting Li abundances given in Table 10 
are estimated to be affected by errors of the order of 0.4-0.6 
dex. As an example, the changes in the derived Li abundance 
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Fig. 8. Best model fit and observed spectrum in the region 7670-7730 A (around the K I line at 7699 A) for the star IRAS 
1 1081-4203. The synthetic spectrum obtained for solar K abundance (log e(K) = 5.12 dex) is shown. The parameters of the best 
model atmosphere fit are indicated in the top label. Note the radial velocity shift, which is interpreted as a result of the partially 
circumstellar origin of the K I line. The narrow features seen at ~7675-77 A, which are not adjusted by the model, are telluric 
lines. 



induced by slight variations of each of the atmospheric param- 
eters used in our modelling for IRAS 15255+1944 are shown 
in Table 1 1 . 

Note, however, that the estimated errors do not reflect pos- 
sible non-LTE efifects, dynamics of the stellar atmosphere, or 
errors in the model atmospheres or in the molecular/atomic 
linelists themselves. In particular, over-ionization and over- 
excitation of Li is predicted to occur under non-LTE conditions 
both in C-rich and O-rich AGB stars (Kiselman & Plez 1995; 
Pavlenko 1996) although the effect is expected to be more pro- 
nounced in metal-deficient stars (i.e. MC AGBs). Use of LTE 
in Li-rich AGB stars is likely to result in underestimation of 
the Li abundance (Kiselman & Plez 1995; Abia, Pavlenko & 
Laverny 1999). 

Concerning s-process elements, the non-detection of the 
ZrO molecular bands at 6474 A and 6495 A is found to be con- 
sistent with upper limits for the Zr abundance around [Zr/Fe] = 
0.25 ± 0.25 dex with respect to the solar value of log e(Zr) = 2.6 
dex. The synthetic spectra predict detectable ZrO bandheads at 
6474 A and 6495 A for Zr abundances above 0.00-0.25 dex 
for Teff > 3000 K and above 0.25-0.50 dex for T^b < 3000 K. 
Actually, the fits are also reasonably good even if we do not 
include any ZrO in the synthesis. We found the same results 
both in the Li non-detected stars and in the Li detected ones. 
This result is in strong contrast with the higher Zr abundances 
([Zr/Fe] > 0.5) found in Galactic MS, S-stars and in massive O- 
rich AGB stars in the MCs. The effect of varying the Zr abun- 
dance for the Li detected star IRAS 11081-4203 is displayed 
in Figure 11, as an example. The best model spectrum fit for 



the Galactic S-star IRAS 10436-3454 is also shown in Figure 
12 for comparison. As we can see, strong ZrO bandheads are 
detected in IRAS 10436-3454, while these are completely ab- 
sent from the spectrum of IRAS 11081-4203. According to 
our models, a modest Zr enhancement (with respect to the solar 
value) would be enough for the ZrO molecular bands to show 
up in our sample stars despite the strong veiling produced by 
the TiO molecule (see Figure 11). An even lower Zr abundance 
is enough to produce strong ZrO bands in S-stars, where the 
C/O ratio is ~1 because of the weaker TiO veiling, as we can 
see in the spectrum of IRAS 10436-3454 displayed in Figure 
12. This also favours the visibility of the Zr 1 atomic lines in the 
7400-7600 and 8 100-8 150 A spectral regions in these stars, as 
has already been mentioned. 

5. Discussion 

5.1. Absolute luminosities 

In constrast to the studies made in the past on AGB stars in 
the Magellanic Clouds, for which a common distance can be 
assumed leading to relatively well determined absolute lumi- 
nosities, the estimation of absolute luminosities for the stars 
in our Galactic sample is very difficult, mainly because of the 
large uncertainties involved in the determination of distances 
within our Galaxy.^ As a first approach, we can try to esti- 

' Note that, in addition, AGB stars are also strongly variable. As an 
example, Engels et al. (1983) found an M(,oi variation of 2 magnitudes 
between the minimum and maximum light (from -5.0 to -7.1 mag 
respectively) in the Galactic GH/IR star OH 32.8-0.3. 
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Fig. 9. Best model fit and observed spectrum in the region 7775-7835A (around the Rb I line at 7800 A) for the star IRAS 
1 108 1 -4203. The synthetic spectra obtained for solar Rb abundance (log e(Rb) = 2.6 dex) and shifted +0.5 dex (log e(Rb) = 3.1 
dex) are shown. The parameters of the best model atmosphere fit are indicated in the top label. Again, the radial velocity shift 
observed suggests a partially circumstellar origin of the Rb 1 line. 



mate luminosities using the period-luminosity relationship (for 
those sources with a well determined period). If the P-L rela- 
tionship for Galactic Mira variables found by Groenewegen & 
Whitelock (1996) is applied to our sample of OH/IR stars, we 
obtain Mboi around — 6 for those stars with periods of 400- 
600 days while an Mboi of — 5.5 is derived for periods around 
~300 days. Similar values were found among the more massive 
AGB stars in the MCs. However, if we take those stars with the 
longer periods in the sample (> 1000 days), we derive Mboi 
< -7, which seems unrealistically high. This suggests that the 
period-luminosity relationship may not hold for the most ex- 
treme OH/IR stars (Wood, Habing & McGregor 1998). On the 
other hand, there is recent observational evidence for the exis- 
tence of low metallicity (e.g. in the LMC) HBB AGB stars with 
luminosities brighter than the predictions of the core-mass lu- 
minosity relation which have been atributed to an excess flux 
from HBB (Whitelock et al. 2003). A similar effect could ex- 
plain the detection of AGB stars with such high Mboi also in 
our Galaxy. 

5.2. Progenitor masses 

Neither is determining the progenitor masses of our Galactic 
O-rich AGB stars a simple task. Several observational param- 
eters, however, such as the OH maser expansion velocity or 
the variability period, have been proposed in the past as useful 
distance-independent mass indicators for this class of stars, and 
we will make use of them in the following. 

The assumption is based on the fact that the group of OH/IR 
stars with longer periods and larger expansion velocities shows 



a Galactic distribution that corresponds to a more massive pop- 
ulation (Baud et al. 1981; Baud & Habing 1983; Chen et al. 
2001; Jimenez-Esteban et al. 2005b) suggesting that the peri- 
ods of variability and the expansion velocities of the circum- 
stellar envelopes of AGB stars must be closely correlated with 
the progenitor mass. 

Actually, stars in our sample follow the expected trend, as 
is illustrated in Figure 13. Although the limited number of stars 
considered prevents us from deriving any statistical conclusion 
based on our data alone, it seems clear that the galactic lati- 
tude distribution of the stars in our sample becomes narrower 
as a function of Vexp(OH), as one would expect if the higher 
expansion velocities do correspond to a more massive popula- 
tion. A similar trend is observed when the variability period is 
analysed (not shown). 

Figure 14 (left panel) shows the distribution of the sources 
in our sample on the IRAS two-colour diagram [12] - [25] vs 
[25] - [60] for which OH expansion velocities are available. 
As we can see in this diagram, on average, the sources with in- 
termediate OH expansion velocities, between 6 and 12 km s~\ 
show redder [12] - [25] and [25] - [60] colours than those with 
low OH expansion velocities (i.e. below 6 km s '), which im- 
plies thicker circumstellar envelopes or larger mass loss rates, 
but they are bluer than the subsample of stars with OH expan- 
sion velocities beyond 12 km s Redder IRAS colours can 
be interpreted as the consequence of a more advanced stage 
on the AGB if we assume that mass loss rate increases along 
this evolutionary phase. Alternatively, they may also be inter- 
preted as corresponding to a more massive population of AGB 
stars since in principle they should be able to develop thicker 
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Fig. 10. High resolution optical spectra around the 7565 A region of two stars in our sample where we show the lack of s-process 
element atomic lines (e.g. from Zr and Nd) in comparison with the Galactic S-star IRAS 10436-3459. IRAS 05027-2158 and 
IRAS 19361-1658 show lithium lines with very different strength. The positions of some atomic lines of Zr, Nd and Fe are 
indicated. It seems that only the Fe I line at ~7586 A is clearly detected in all stars. 



circumstellar envelopes even at a relatively early stage in their 
AGB evolution. A similar behaviour is observed when the vari- 
ability period is analysed (Figure 14; right panel), although the 
number of stars with known periods is significantly smaller. In 
this case we can see that the reddest sources in the diagram be- 
long to the subgroup showing the longer periods, in excess of 
700 days. 

The distribution of variability periods and OH expansion 
velocities is found to be remarkably different for each of the 
three subgroups identified in our sample. Histograms showing 
the results obtained are shown in Figures 15 and 16, respec- 
tively. The same results are presented in a different way in Table 
12, where the distribution of sources among the different sub- 
types is presented as a function of their OH expansion velocity 
and variability properties. 

In Table 12 we can see that most of the AGB stars with 
the lower OH expansion velocities (vexp(OH) < 6 km s~^) are 
among those undetected in the Li I line (92%) while 44% of 
stars with OH expansion velocities between 6 and 12 km s"' 
(but still observable in the optical) are detected as Li-strong 
sources. However, there are also stars with non-detection of 
lithium (24%) and stars too red or not found (32%) within this 
group. Finally, a majority of the stars in the group of heavily 
obscured sources (too red or not found) shows OH expansion 
velocities higher than 12 km s~' (66%). 

Similar statistics are obtained for the variabiUty period. 
Sources with relatively short periods, i.e. below 400 days, are 
predominantly non-Li detections (84%), while almost half of 
the stars showing periods between 400 and 700 days are Li 
detected (46%). Again, the same observational problems arise 



when we try to analyse the stars with the longer periods {P > 
700 days), since most of them are completely obscured in the 
optical (88%). 

5.3. Period and OH expansion velocity versus HBB 

In Figure 17 we can see that there is actually a clear corre- 
lation between the Li abundance and the variability period in 
our target stars. Remarkably, for periods below 400 days (log 
P < 2.6), only one star out of 12 is Li detected. In contrast, al- 
most all stars with P > 500 days (log P > 2.7) are Li-rich, the 
only exception being IRAS 18050-2213 (= VX Sgr), which is 
a peculiar, very cool semi-regular pulsating star, not showing a 
proper Mira variability in the last decade, which some authors 
identify as a massive red supergiant, and not as an AGB star 
(Kamohara et al. 2005). 

A similar correlation, although not so clear, is found when 
the expansion velocities derived from the OH maser measure- 
ments are considered (see Figure 18). In this case we can see 
that the maximum Li abundance increases with the OH expan- 
sion velocity, suggesting that the most massive AGB stars can 
experience a higher lithium enrichment. However, the more in- 
teresting result is the non-detection of Li in any star of the sam- 
ple with Vexp(OH) < 6 km s-^ IRAS 18050-2213 is again one 
of the two outliers showing no Li and a very high expansion 
velocity, while the other one is IRAS 13442-6109, for still un- 
known reasons. 

At present, it is generally accepted that Uthium production 
in massive O-rich AGB stars is due to the activation of hot bot- 
tom burning at the bottom of the convective mantle, which pre- 
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Table 10. Spectroscopic Teff and Li abundances' derived for the subgroup of Li detected (left) and Li non-detected stars (right) 
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' The uncertainty in the derived Li abundances is estimated to be around ±0.4-0.6 dex 
^ Possible double lined spectroscopic binary 
' The lithium line has a P cygni-type profile 

Possible non-AGB star 
' Observed Li abundances in two different runs 

* The abundance value must be treated with some caution because the line is resolved in two components (circumstellar and stellar) and the 
abundance estimate corresponds to the photospheric abundance needed to fit the stellar component 
** The S/N ratio at 6708 A is very low to derive any reUable Teff and/or Li abundance estimate 



vents the formation of luminous C-rich AGB stars (see details 
in Section 1). If the detection of lithium is taken as a signature 
of HBB and the period and the OH expansion velocity are ac- 
cepted as valid mass indicators, our results indicate that there 
are no HBB stars among the subgroup of AGB stars in our sam- 
ple with periods lower than 400 days and Vexp(OH) < 6 km s \ 
while most of the AGB stars with periods beyond 500 days and 
high expansion velocities are quite probably developing HBB. 

5.4. Lithium production and HBB model predictions 

Modelling of HBB started thirty years ago with envelope mod- 
els including non-instantaneous mixing coupled to the nuclear 



evolution (Sackmann, Smith & Despain 1974). More recent 
models (Sackmann & Boothroyd 1992; Mazzitelli, D'Antona 
& Ventura 1999, hereafter SB92 and MDV99, respectively) 
seem to work well in reproducing the lithium production ob- 
served in massive O-rich AGB stars in the MCs. Although the 
precise determination of the lithium production depends on the 
input physics (stellar mass, metallicity, mass loss rate, over- 
shooting, etc.), the strongest dependence is with the convection 
model assumed. SB92 models predict lithium production in in- 
termediate mass stars assuming the so-called "mixing length 
theory" (MLT) framework for convection. But these models re- 
quire a fine tuning of the mixing length parameter a. MDV99 
explored lithium production in intermediate mass stars accord- 
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Table 11. Sensitivity of the derived Li abundances (in dex) to slight changes in the model atmosphere parameters for IRAS 
15255+1944 
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Fig. 11. Synthetic and observed spectra in the region 6460-6499 A for the star IRAS 1 1081-4203, one of the Li detected stars 
in our sample. The synthetic spectra corresponding to [Zr/Fe] = +0.0, +0.5, and +1.0 dex (or log e(Zr) - 2.6, 3.1, and 3.6 dex, 
respectively) are shown. The parameters of the best model atmosphere fit are indicated in the top label. 

Table 12. Distribution of sources among the different subgroups identified in the sample as a function of their OH expansion 
velocity and variability properties. 



Period (days) Observed Li undetected Li detected Too red/not found 

<400 12 10 (84%) 1 (8%) 1 (8%) 

400-700 26 4(15%) 12(46%) 10(39%) 

>700 17 1 (6%) 1 (6%) 15 (88%) 



Vexp(0/f) (km s-') 


Observed 


Li undetected 


Li detected 


Too red/not found 


<6 


13 


12 (92%) 


(0%) 


1 (8%) 


6-12 


25 


6 (24%) 


1 1 (44%) 


8 (32%) 


>12 


53 


10 (19%) 


8 (15%) 


35 (66%) 



ing to a more modern treatment of turbulent convection, known 
as the "Full Spectrum of Turbulence" (FST) model (D'Antona 
& Mazzitelli 1996, and references therein). Actually, MDV99 
models are able to reproduce the strong Li enhancements ob- 
served in massive O-rich LMC AGB stars in a satisfactory way 
and provide a theoretical explanation to the existence of a few 
sources in the LMC with extreme luminosities of up to Mboi = 



-7.3 and -7.6, which are known to be long period, obscured 
AGB stars (Ventura, D'Antona & MazzitelU 2000). However, 
these models have not yet been tested against their Galactic 
analogues. 

For solar metallicity, MDV99 predict lithium production as 
a consequence of the activation of HBB for stellar masses > 
4.5 Mq without core overshooting and for M > 4 M© including 
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T,„=3aOQ K, logg-0.35, z-0.0, C/O-0.9, (^3 km s-', FWHM^400mA 
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Fig. 12. Synthetic and observed spectra in the region 6460-6499 A for the Galactic S-star IRAS 10436-3454. The synthetic 
spectra corresponding to [Zr/Fe] - +0.0, +0.5 and +1.0 dex (or log e(Zr) = 2.6, 3.1, and 3.6 dex, respectively) are shown. The 
parameters of the best model atmosphere fit are indicated in the top label. Note the strong overabundance of [Zr/Fe] = +1.0 dex 
(log e(Zr) - 3.6 dex) needed to fit the observations. 
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Fig. 13. Galactic latitude distribution of the sources in our sample as a function of their OH expansion velocities: Vexp(OH) < 6 
km s ' (left panel); 6 < Vexp(OH) < 12 km s ' (middle panel); and Vexp(OH) > 12 km s ' (right panel). 



core overshooting. The amount of lithium produced during the 
AGB does not depend on the assumed initial lithium abundance 
as the star structure loses all memory of the previous history of 
lithium at the beginning of the AGB phase. As an example, ac- 
cording to this model the lithium surface abundance drops to 
log e(Li) ~ 0.3 for a 6 Mq star after the second dredge-up, at 
the beginning of the AGB phase. If overshooting from below 
the convective envelope is allowed, an important decrease in 
the number of thermal pulses is predicted by the model and 
lithium production takes place before the first thermal pulse. 
The variation with time of the lithium abundance, total lumi- 
nosity and temperature at the base of the convective envelope 
(^bce) for different masses is shown in Figure 15 of MDV99. A 
stronger Li-overabundance, faster increase of the luminosity at 



the beginning of the AGB phase and higher Tbce are predicted 
when the mass of the star is increased. 

Mass loss rates and metallicities play also an important role 
according to MDV99. The minimum mass needed for the ig- 
nition of HBB decreases with decreasing mass loss rates. In 
addition, a low mass loss rate leads to a longer run of ther- 
mal pulses, and to a lithium abundance which remains large 
for the entire thermal pulsing AGB (TP-AGB, hereafter) life- 
time. In contrast, for high mass loss rates, the HBB process is 
stopped only after a few pulses (~5 for a 6 Mq star). Models 
with strong mass loss show oscillations of the Li abundance by 
orders of magnitude on short timescales of ~10'* years (even 
within a single thermal pulse!). The temperature at the bot- 
tom of the convective envelope, Tbce, varies so much that the 
hthium production is temporarily stopped while the lithium al- 
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Fig. 14. Distribution of the sources included in our sample in the IRAS two-colour diagram [12] - [25] vs [25] - [60] showing 
the location of the stars in our sample as a function of the OH expansion velocities (left panel) and variability periods (right 
panel). The continuum line is the "O-rich AGB sequence" (see text), which indicates the sequence of colours expected for O-rich 
AGB stars surrounded by envelopes with increasing thickness and/or mass loss rates (Bedijn 1987). 
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Fig. 15. Distribution of variability periods for the Li non-detected (left panel), Li detected (middle panel) and too red/not found 
(right panel) sources. 
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Fig, 16. Distribution of OH expansion velocities for the Li undetected (left panel), Li detected (middle panel) and too red/not 
found (right panel) sources. 
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Fig. 17. Observed Li abundance vs variability period. Upper 
limits to the Li abundance are shown with black triangles and 
marked with arrows. Abundance values derived from Li I lines 
which are resolved in two components (circumstellar and stel- 
lar) are shown with red triangles and correspond to the pho- 
tospheric abundance needed to fit the stellar component. The 
dashed vertical line corresponds to a variability period of 400 
days (see text). 
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Fig. 18. Observed Li abundance vs OH expansion velocity. 
Upper limits to the Li abundance are shown with black trian- 
gles and marked with arrows. The abundance values derived 
from Li I lines which are resolved into two components (cir- 
cumstellar and stellar) are shown with red triangles and corre- 
spond to the photospheric abundance needed to fit the stellar 
component. The dashed vertical line corresponds to an OH ex- 
pansion velocity of 6 km s"' (see text). 



ready produced has already been diluted by convection. When 
Tbce increases, the lithium abundance increases again. The con- 
sequence is that, at least for 20% of the time, there is negligible 
lithium in the envelope (at a level of approximately log e(Li) < 
1.0). This implies that the distribution of Li abundances derived 
from the observations can only be analysed in a statistical way. 

On the other hand, the activation of HBB takes place at a 
lower mass limit of only 3.0-3.7 Mq (depending on the mass 
loss rate considered) at the metallicity of the LMC (Ventura, 
D'Antona & Mazzitelli 2000). Low metallicity models predict 
larger values of Tbce and higher luminosities, and as a conse- 
quence, higher lithium production for a given stellar mass. In 
addition, this accelerates Li-production and a maximum value 
of log e(Li) ~ 4.3 is reached in a shorter timescale with respect 
to the solar metallicity case. The variation with time of the sur- 
face lithium abundance, stellar luminosity and T^ce, computed 
for a 6 Mq model and different metallicities is shown in Figure 
10 of MDV99. Lower metallicity models esentially produce 
lithium faster due to the larger temperature at the base of the 
convective envelope, but show high Li abundances for a shorter 
time. 

5.5. Theory versus observations 

Globally considered, our results indicate that below a certain 
value of the period and the OH expansion velocity our stars do 
not show any Li enhancement. However, not all the stars above 
these values show detectable Li. The strong IR excesses asso- 
ciated with most of the stars in our sample suggest that they are 
actually at the end of their evolution as AGB stars. Surprisingly, 
none of the stars in our sample shows a strong s-process en- 
hancement (e.g. from Zr, Y, Sr) as would have been expected if 
they were stars that have experienced a significant third dredge- 
up during their AGB evolution. In the case of the more massive 
stars in our sample one could argue that the mass loss is so 
strong at the high metallicity of our Galaxy (compared to the 
MC metallicity) that their AGB lifetime may become shortened 
significantly, so that they will end this phase experiencing too 
few thermal pulses before they completely lose their envelope. 
However, the lack of s-process enhancement is also observed 
in the group of stars in which no lithium has been detected. 

In order to interpret these results, several scenarios can be 
proposed: 

i) all the stars in the sample are HBB stars, and thus, consid- 
erably massive (M > 4 Mq). The difference between the 
Li detected and the Li undetected ones would be just that 
the undetected ones are in a very early stage as AGB stars. 
This hypothesis would explain the absence of s-process en- 
hancement (e.g. Zr, Y, Sr) in all the stars in the sample, 
but it would not be consistent with the non-detection of 
lithium in the group of stars with the shorter periods and 
lower OH expansion velocities. HBB models predict a con- 
siderable Li enhancement just after the arrival at the AGB. 
Statistically, there should be no difference in the behaviour 
of the Li abundances for any star in the sample, and this is 
not observed. 
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the sample is composed by a combination of low mass 
(M < 1.5 Mq) and high mass (M > 4 Mq) stars. In this sce- 
nario, the optically bright O-rich AGB stars in our sample 
with the shorter periods and lower OH expansion velocities, 
undetected in Li, would belong to the population of low 
mass stars in the Galaxy that do not experience HBB and 
evolve all the way along the AGB as O-rich stars. However, 
these stars are expected to evolve very slowly and develop 
a large number of thermal pulses during the AGB. At the 
end of their AGB evolution they are expected to become 
strongly enhanced in s-process elements as a consequence 
of the dredge-up of processed material induced by the many 
thermal pulses experienced. But the Li undetected stars in 
our sample do not show this enhancement in s-process ele- 
ments. The only way to reconcile this result with their sta- 
tus as low mass stars would imply assuming that they are 
at a very early stage as AGB stars, which again seems to be 
incompatible with the strong IR excess detected by IRAS. 
Such a strong excess is only expected in low mass stars at 
the very end of the AGB. Alternatively, they could experi- 
ence less thermal pulses than predicted by the models or the 
dredge-up process could be very inefficient. But this is not 
what we observe in other low and intermediate mass stars 
in the Galaxy, such as the S-stars and many C-rich stars, 
which show the strong s-process enhancement that is not 
observed in our stars. Therefore, we must conclude that the 
hypothesis of these Li undetected stars being identified as 
low mass AGB stars seems to be implausible, 
all the stars in the sample are intermediate mass stars (M > 
3 Mq) but only a fraction of them, the more massive ones 
(M > 4 Mq), experience HBB. This scenario would ex- 
plain why the sources with lower periods and OH expan- 
sion velocities in our sample are undetected in the Li I line, 
while at the same time, they do not show any s-process en- 
hancement. These stars, although slightly below the mini- 
mum mass needed to develop HBB, would also experience 
a very strong mass loss rate and only a few thermal pulses 
(~5-10) before the early termination of the TP-AGB as a 
consequence of the strong mass loss. The initial envelope 
of these stars would also be quite massive and thus a slow 
change of the C/O ratio is expected because of the huge di- 
lution of the dredged up material. But how can we explain 
the large spread of Li abundances observed in the group of 
sources in our sample with longer periods and larger OH 
expansion velocities? These are all supposed to be HBB 
stars. And HBB models predict a large fluctuation in the Li 
abundance as the star evolves along the TP-AGB (see fig- 
ure 15 of MDV99 again) between two consecutive pulses. 
Actually, the timescale of the Li production phase in the 
TP-AGB is predicted to be of the order or slightly smaller 
than the timescale between pulses (~10'* years for a mas- 
sive star with M > 4-Mq). In addition, as the AGB evo- 
lution proceeds, the surface Li abundances decrease again 
very rapidly as a consequence of the ^He becoming almost 
completely burned in the envelope and thus the Li produc- 
tion ceases (e.g. MDV99; Forestini & Charbonnel 1997). 
This would be consistent with the detection of strong Li 
abundances in ~50% of the stars in our sample for which 



we suspect that HBB is active. In this scenario, the strongly 
obscured sources for which we were not able to obtain an 
optical spectrum could also be interpreted as a population 
of high mass stars developing HBB with even more extreme 
observational properties. They are sources with longer pe- 
riods (> 700 days) and larger OH expansion velocities (> 
12 km s ') which may have developed optically thick cir- 
cumstellar envelopes as a consequence of the huge mass 
loss experienced in the AGB. These stars must also be de- 
tectable as Li-rich at least for some time during their AGB 
evolution but, unfortunately, they are so strongly reddened 
in the optical range that their study is only possible in the 
infrared, where no strong Li lines are present. We propose 
that the HBB nature of these obscured stars can be eluci- 
dated by determining their '^C/'^C ratios, which can also 
be used as an HBB indicator, through near infrared spec- 
troscopy. 

5.6. Comparison with tine 0-ricin AGB stars in tlie 
Mageiianic Ciouds 

The situation is quite dififerent for the AGB stars in the MCs. 
In contrast with their Galactic analogues, the more luminous 
AGB stars in the MCs (with -7 < Mboi < -6) are O-rich stars 
showing s-process element enhancements. In addition, a higher 
proportion of them (~80% compared to ~50% in our Galaxy) 
also shows lithium overabundances. The enrichment in lithium 
indicates that they are also HBB stars, but why are these stars 
also enriched in s-process elements? 

The answer to this question must be related to the difl'erent 
metallicity of the stars in the MCs with respect to our Galaxy. 
Actually, theoretical models predict a higher efficiency of the 
dredge-up in low metallicity atmospheres (e.g. Herwig 2004) 
with respect to those with solar metallicity (e.g. Lugaro et al. 
2003). 

Moreover, there is observational evidence that lower metal- 
licity environments are also less favourable to dust production, 
as is suggested by the smaller number of heavily obscured AGB 
stars in the MCs (Trams et al. 1999; Groenewegen et al. 2000), 
compared to our Galaxy. This is supported by the lower dust- 
to-gas ratios derived by van Loon (2000) in the few obscured 
AGB stars in the MCs for which this analysis has been made. 

If mass loss is driven by radiation pressure on the dust 
grains, it might be less efficient with decreasing metallicity 
(Willson 2000). In that case, longer AGB lifetimes would 
be expected, which could increase the chance of nuclear- 
processed material reaching the stellar surface. This would also 
explain the larger proportion of luminous C-rich stars (up to 
Mboi ~ -6) observed in the MCs (Plez, Smith & Lambert 1993; 
Smifli & Lambert 1989, 1990a; Smifli et al. 1995) with respect 
to the Galaxy. The slow evolution predicted for AGB stars in 
the MCs as a consequence of the less efficient mass loss leaves 
time for more thermal pulses to occur during the AGB life- 
time and, therefore, a more effective dredge-up of s-process el- 
ements to the surface can be expected before the envelope is 
completely gone at the end of the AGB. This would explain 
why even the more massive stars in the MCs show a strong s- 
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process enrichment in contrast to their Galactic counterparts. 
In our Galaxy the only AGB stars showing a similar overabun- 
dance in s-process elements seem to be the result of the evo- 
lution of low to intermediate mass stars (M < 1.5-2.0 Mg), 
while no, or very little, s-process enhancement is observed in 
Galactic AGB stars with higher main sequence masses. 

In addition, the lower critical mass (which also decreases 
with decreasing mass loss rate, see Section 5.4) needed to de- 
velop HBB (M > 3 Mq at the metallicity of the LMC, com- 
pared to the ~4 Mq limit in our Galaxy) favours the simulta- 
neous detection of s-process elements and Li enrichment in a 
larger number of AGB stars in the MCs. In our Galaxy, the Li- 
rich sample of AGB stars is restricted to the fraction of stars 
with main sequence masses M > 4 Mq. In contrast to their 
MC counterparts, these stars evolve so rapidly (because of the 
strong mass loss) that there is no time for a significant enhance- 
ment in s-process elements, as the results here presented seems 
to demonstrate. 

Obviously, more refined s-process calculations taking into 
account the efi'ects of HBB and mass loss, using a large grid of 
stellar evolutionary models with different masses and metallic- 
ities are strongly encouraged. 

5.7. Evolutionary connections 

The well known existence of at least two different chemical 

branches of AGB stars (C-rich and O-rich) and the dependence 
of this chemical segregation on the main sequence mass and 
on the environmental conditions in which these stars evolve 
(such as the metallicity) is a long debated issue in stellar evolu- 
tion. Moreover, the connection with the chemical composition 
observed in Galactic post- AGB stars and/or well evolved PNe 
(which are expected to be the result of the evolution of AGB 
stars with a wide variety of masses) has never been addressed 
in detail from an overall perspective. 

Stars classified as post- AGB in the literature are mainly low 
mass stars, according to their wide Galactic distribution. They 
are thought to represent the fraction of stars which evolve so 
slowly that they are still visible for some time as stars with 
intermediate spectral types before they become PNe. Many of 
these stars are C-rich and show strong s-process element en- 
hancements [s-process/Fe] ~ +1.5 (van Winckel & Reyniers 
2000; Reddy, Baker & Hrivnak 1999; Klochkova et al. 1999), 
in agreement with what we would expect if they were the fi- 
nal product of the evolution of low mass AGB stars. However, 
other post-AGB stars do not show the characteristic signa- 
tures of the third dredge-up. In the non-enriched sources the 
s-process element abundances can be as low as [Zr/Fe] < 0.2 
(Luck, Bond & Lambert 1990; van Winckel 1997). These ob- 
jects might be stars with even lower masses (M ~ 0.8-1.0 Mq) 
which have evolved off the AGB before the third dredge-up 
could occur (Marigo, Girardi & Bressan 1999). Some of these 
sources are high galactic latitude, hot B-type post-AGB stars 
(e.g. Conlon et al. 1993; Mooney et al. 2001) that show strong 
carbon deficiencies. In the more extreme cases they will prob- 
ably never become PNe. 



On the other hand, and as we have shown in this article, 
the more massive AGB stars may also show a similar lack of s- 
process elements, but in this case as a consequence of the rapid 
evolution induced by the strong mass loss. These stars are more 
difficult to observe in the optical as they are expected to evolve 
as heavily obscured sources from the TP- AGB to the PN stage 
but may be the progenitors of most Galactic PNe. 

Actually, a similar chemical segregation is also observed in 
Galactic PNe. Classically, they are classified as type I, II or III 
as a function of their chemical abundance properties (Peimbert 
1978). They are expected to cover a wide range of progenitor 
masses and thus show different chemical signatures as a conse- 
quence of their previous passage along the AGB. 

In particular, the higher mass fraction of stars (with M > 
4 Mq) identified in this paper as stars developing HBB and 
showing a strong lithium enhancement are expected to arrive at 
the PN stage still as heavily obscured stars. These stars would 
evolve into N-rich type I PNe, which are characterized by their 
large He abundances (He/H > 0.10) and N/O ratios (Manchado 
2003, 2004; Mampaso 2004, and references therein). They are 
thought to represent the fraction of the more massive PNe in 
the Galaxy, as is confirmed by their strong concentration at 
low galactic latitudes. Their strong N overabundances would 
be consistent with their identification as PNe resulting from the 
evolution of HBB AGB stars. 

6. Conclusions 

As a result of a high resolution spectroscopic survey carried 
out in the optical for a large sample (102) of massive Galactic 
O-rich AGB stars, we have detected the Li I resonance line at 
6708 A in 25 stars, while 32 stars did not display this Li I line. 
The rest of the stars observed (45 stars) were heavily obscured 
by their thick circumstellar envelopes and, as such, they were 
too red or not found at optical wavelengths. 

We have classified the sample into three subgroups on the 
basis of their OH expansion velocity and variability period. 
From their relative distribution in the IRAS two-colour dia- 
gram [12] - [25] vs [25] - [60] and Galactic distribution, we 
conclude that they must represent populations of Galactic O- 
rich AGB stars with different progenitor masses. 

Combining MARCS model atmospheres and synthetic 
spectroscopy with extensive line lists we have been able to de- 
rive the fundamental stellar parameters (Teff , log g, C/O, ^, etc.) 
and obtain the Li and Zr abundances in those stars for which an 
optical spectrum was obtained. Our chemical abundance analy- 
sis shows that half of stars show Li overabundances in the range 
between log e(Li) ~ h-0.5 and +3.0 dex. This is interpreted as 
a signature of the activation of HBB, confirming that they are 
actually massive AGB stars (M > 4 Mq). However, these stars 
do not show any zirconium enhancement (taken as representa- 
tive for the s-process element enrichment), indicating that they 
behave differently from the more massive (and luminous) AGB 
stars in the MCs. 

Assuming that the variabiUty period and the OH expansion 
velocity can be taken as distance-independent mass indicators 
and by comparing our results with the theoretical predictions of 
HBB and nucleosynthesis models we conclude that the O-rich 
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AGB stars in our sample are all considerably massive (M > 3 
Mq) even those non-detected in lithium, but only the more mas- 
sive ones (M > 4 Mq; with periods longer than 400 days and 
Vexp(OH) > 6 km s experience HBB. The lack of lithium en- 
richment found in some of the sources thought to belong to the 
group of more massive AGB stars in the sample is explained 
as the consequence of the timescale of the Li production phase 
being of the order of or sUghtly smaller than the interpulse time 
(~10"* years) or, alternatively, by the ^He exhaustion in the en- 
velope interrupting the Li production, in agreement with the 
predictions made by the models. 

The group of stars displaying the most extreme observa- 
tional properties (with periods sometimes longer than 700 days 
and Vexp(OH) > 12 km s"^) are beUeved to represent the more 
massive AGB stars in our Galaxy. Unfortunately, these stars 
are strongly obscured by their thick circumstellar envelopes 
and we could not carry out any type of analysis in the optical 
range. This very high mass population should also experience 
HBB and show strong Li enhancement but no s-process ele- 
ment overabundances. We propose that the HBB status of these 
obscured stars can be determined by measuring their '^C/'^C 
ratios in the near-infrared, which can also be used as a HBB 
indicator. 

Our results suggest that the dramatically different abun- 
dance pattern found in AGB stars belonging to the MCs and 
to our Galaxy can be explained in terms of the different metal- 
licity conditions under which these stars evolved. This consti- 
tutes strong observational evidence that the chemical evolution 
of massive AGB stars is strongly modulated by metallicity. 

Globally considered, our results can also be used to es- 
tablish evolutionary connections between TP- AGB stars, post- 
AGB stars and PNe. In particular, we suggest that the HBB 
AGB stars identified in our sample must be the precursors of 
N-rich type I PNe. 
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